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Abstract 

We present a complete calculation of the higher-order perturbative QCD corrections to in- 
elastic photoproduction of J particles. A comprehensive analysis of total cross sections and 
differential distributions for the energy range of the fixed-target experiments and for inelastic 
J /ill photoproduction at HERA is performed. The cross section and the J/i/' energy spectrum 
are compared with the available photoproduction data including first results from HERA. This 
analysis will not only provide information on the gluon distribution of the proton but appears 
to be a clean test for the underlying picture of quarkoniuni production as developed so far in 
the perturbative QCD sector. 
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1 Introduction 



The measurement of the gluon distribution in the nucleon is one of the important goals of lepton- 
nucleon scattering experiments. The classical methods exploit the evolution of the nucleon struc- 
ture functions with the momentum transfer and the size of the longitudinal structure function. 
With rising energies, however, jet physics and the production of heavy quark states become im- 
portant complementary tools. Heavy flavour production in lepton- nucleon scattering is dominated 
by photon-gluon fusion and can thus yield direct information on the gluon distribution in the nu- 
cleon G{x, Q^). Besides open charm and bottom production, the formation of J/ip bound states in 
inelastic photoproduction experiments 

j + M ^ J/i; + X (1) 

provides an experimentally attractive method since J/ip particles are easy to tag in the leptonic 
decay modes. 

The production of heavy quarks in high energy photon-proton collisions can be calculated in 
perturbative QCD. The mass of the heavy quark, uiq ^> Aqcd, acts as a cutoff and sets the scale 
for the perturbative calculations Q. However, the subsequent transition from the colour-octet QQ 
pair to a physical quarkonium state introduces non-perturbative aspects. Two different mechanisms 
of bound state formation have been employed in previous analyses (for a recent review see Ref. [Q]): 

(i) The local duality approach Q assumes that the colour-octet QQ pair rearranges itself into 
a colour-singlet bound state by the emission of non-perturbative soft gluons. According to the 
arguments of semi-local duality, one averages over all possible quarkonia states by integrating the 
perturbative cross section for inclusive QQ production from the quark threshold (= 2mQ) to the 
physical threshold for the production of a pair of heavy- light mesons (= 2m d for the cc system) . The 
probability to generate a particular state, e.g. J/ip, depends on different dynamical details of the 
production mechanism and cannot be absolutely predicted in this model. Another serious drawback 
of the duality approach is the fact that higher-order QCD corrections cannot be included since there 
is no unique way to decide what part of the radiatively emitted gluons are to be considered as a part 
of the bound state. Although dual models might describe some qualitative features of quarkonium 
production, they do not allow for a quantitative prediction and will therefore not be discussed in 
the present context. 

(ii) In the colour- singlet ( CS) mechanism the quarkonium state is described by a colour- 
singlet QQ pair with the appropriate spin, angular-momentum and charge conjugation quantum 
numbers. In the static approximation, in which the motion of the charm quarks in the bound state 
is neglected, the production cross section factorizes into a short distance matrix element which 
describes the production of a QQ pair within a region of size l/mg, and a long distance factor that 
contains all the nonperturbative dynamics of the bound state formation. The short distance cross 
section can be calculated as a perturbative expansion in powers of the strong coupling constant 
as^mq), evaluated at a scale set approximately by the heavy quark mass, while the long-distance 
factor is related to the leptonic width. 

A rigorous framework for treating quarkonium production and decays has recently been de- 
veloped in Ref.0 (see also Ref.[^). The so-called factorization scheme is based on the use of 
non-relativistic QCD (NRQCD) |^ to separate the short distance parts from the long-distance 
matrix elements. The factorization approach explicitly takes into account the complete structure 
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of the quarkonium Fock space. For the production of S-wave quarkonia, hke J /ip, the colour-octet 
matrix elements associated with higher Fock states like \QQg > are suppressed by a factor of v'^ 
compared to the leading colour-singlet contributions, v being the average velocity of the heavy 
quark in the quarkonium rest frame.§ The NRQCD description of S'-wave quarkonia production 
or annihilation thus reduces to the colour-singlet model in the non-relativistic limit t; — > 0. It has 
been shown in the rigorous analysis of Ref.[^ that the factorization assumption of the CS model is 
correct for any specific /S-wave process in the non-relativistic limit to all orders in as- 

Colour-octet matrix elements can become important if their associated short distance coefficient 
is enhanced compared to the colour-singlet contribution. It has recently been demonstrated that 
large p± charmonium production in hadronic collisions can be accounted for in a satisfactory way 
by including both fragmentation mechanisms as well as higher Fock state contributions |ll,12|. For 



the colour-octet matrix elements associated with the higher Fock states no simple relation exists 
in general between decay and production matrix elements. The corresponding contributions thus 
involve unknown non-perturbative parameters which, in the analyses Refs.|ll|,12], have been fitted 



to the experimental data. Once they have been measured in some production process they can be 
used to predict cross sections for different energies and different beam types. It should however be 
kept in mind that the analyses carried out so far are leading-order analyses. They are therefore 
plagued by large scale dependences and have to rely on the assumption that the perturbative 
expansion of the short-distance coefficients is well behaved and that the higher-order corrections do 
not strongly depend on the specific production mechanism and the collision energy. Moreover, it 
has been argued recently that important higher-twist effects have to be included in the theoretical 



description of charmonium hadroproduction |13]. These effects can yet not be predicted from first 
principles. 

Many channels contribute to the generation of J/ip particles in photoproduction experiments 
[ p!4[ , similarly to the case of hadroproduction. Theoretical interest so far has focussed on two mech- 
anisms for J/tj; photo- and electroproduction, elastic/diffr active [p!^,p!^] and inelastic production 
through photon-gluon- fusion While one expects to shed light on the physical nature of the 
pomeron by the first mechanism, inelastic J/V' production provides information on the distribution 
of gluons in the nucleon The two mechanisms can be separated by measuring the J/tp energy 



spectrum, described by the scaling variable 

z=p-k^/p-ky (2) 

with p, k^^-y being the momenta of the nucleon and J/'0, 7 particles, respectively. In the nucleon 
rest frame, z is the ratio of the J /ip to the 7 energy, z = E^jE^. For elastic/diffractive events z is 
close to one; a clean sample of inelastic events can be obtained in the range z <^ 0.9 [jl^. Reducible 
background mechanisms, such as BB production with subsequent decay into J/v^ particles or 
resolved photon processes might become sizable at HERA energies but can be easily eliminated by 



applying suitable cuts [14|. 



In contrast to the case of hadroproduction, higher Fock state contributions to inelastic J/V' pho- 
toproduction are strongly suppressed compared to the leading colour-singlet mechanism. Colour- 
octet contributions can become important only in the elastic domain, z ~ 1, where the short 
distance coefficient of the colour-singlet amplitude is suppressed by a factor Og. In the inelastic 



^In the case of P-wave quarkonia, colour-singlet and colour-octet mechanisms contribute at the same order in v 
to annihilation rates and production cross sections and must therefore both be included for a consistent calculation 

0- 
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region, z <^ 0.9, a perturbative gluon has to be emitted in the final state for kinematical reasons. 
This gluon carries off colour charge so that a colour-singlet state can be produced without any 
additional short-distance suppression. Since the non-perturbative matrix element associated with 
the colour-singlet state is enhanced compared to the colour-octet matrix element by a factor 
higher Fock state contributions to inelastic J photoproduction can safely be neglected. Fragmen- 



tation mechanisms dominate the production of J/V' particles at large p_L |19| but their contribution 
to the total cross section is very small. As a result, the dominant mechanism for inelastic J/-;/; 
photoproduction is by the colour-singlet state (as assumed in the CS model) and the cross section 
can be predicted without including unknown non-perturbative parameters. 

A comparison of the leading-order predictions with photoproduction data of fixed-target exper- 



iments 1 20, 21 1 reveals that the J/^/; energy dependence AajAz is adequately accounted for in the 
inelastic region z <~ 0.9. The theoretical calculation however underestimates the normalization of 
the measured cross section in general by more than a factor two, depending in detail on the J/'0 



energy and the choice of the parameters |14|. The discrepancy with cross sections extrapolated 
from electroproduction data |2^,23] is even larger. 

The lowest-order approach to inelastic J/V' photoproduction demands several theoretical refine- 
ments: (i) Relativistic corrections due to the motion of the charm quarks in the J/V' bound state; 
(ii) Higher-order perturbative QCD corrections; and last but not least, (iii) Higher-twist effects 
which are not strongly suppressed due to the fairly low charm-quark mass. While the problem of 
higher- twist contributions has not been quantitatively approached so far (see Ref . [|24| for a recent 
discussion), the relativistic corrections have been demonstrated to be under control in the inelastic 
region |]25| , |26|| . Including higher-order QCD corrections is, however, expected to be essential. Their 
contribution in general not only changes the overall normalization of the cross section but can also 
affect the shape of inclusive differential distributions. Expected a 'priori and verified subsequently, 
the QCD corrections dominate the relativistic corrections in the inelastic region, being of the or- 
der of several as(Mj^^) ~ 0.3. In the first step of a systematic expansion, they can therefore be 
determined in the static approach ||^. 

In this work we present the complete calculation of the higher-order perturbative QCD correc- 
tions to inelastic J j'^) photoproduction; first results have already been published in Refs.p7|-p9|]. 
We perform a comprehensive analysis of total cross sections and differential distributions for the 
energy range of the fixed-target experiments and for inelastic J/V' photoproduction at HERA. A 
comparison of the next-to-leading order prediction with the experimental data will not only provide 
information on the gluon distribution of the proton but appears to be a clean test for the underlying 
picture of quarkonium production as developed in the perturbative QCD sector. 

2 The Born cross section 

Inelastic J/V' photoproduction through photon-gluon fusion is described in leading order by the 
subprocess 

7(fci) + 5(^:2) ^JM^)+ 5(^:3) (3) 
as shown in Fig.|l[ Colour conservation and the Landau- Yang theorem |30| require the emission of 



a gluon in the final state. The cross section is generally calculated in the static approximation in 
which the motion of the charm quarks in the bound state is neglected. In this approximation the 
production amplitude factorizes into the short distance amplitude (7 + 5 — > cc + g() , with cc in 



4 



'WWWVWt 











- — ^ 





(a) 



wvwww> 



'WWVWWV 



nrireireisirsir^ 
(b) 



nreirsireireir^ 

(c) 



Figure 1: The leading-order Feynman diagrams contributing to inelastic J/ip photoproduction. 
Additional graphs are obtained by reversing the arrows on the heavy quark lines. 

the colour-singlet state and zero relative velocity of the quarks, and the wave function f{0) of the 
J ftp bound state at the origin: 



M{-f + g^ J/tl^ + g) 



M 



J/V 



ipiO)TT[Us=i,sSPj/^,Mj/^)M{^ + g^cc + g)} . (4) 



The spin projection operator H combines the quark and antiquark spins to the appropriate triplet 

(5) 



states. For negligible binding energy or, equivalently, for nic = one obtains |31] 



ns= 



+ M 



%/2 



where e'^(s^) is the J/V' polarization vector. The coupling strength of the J/V' to the cc-pair is 
specified in terms of the orbital wave function at the origin in momentum space. In leading order, 
the wave function is related to the leptonic width according to 



(6) 



We will describe the cross section in terms of the kinematical (pseudo-) Mandelstam variables 
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where s + t + u = Mj^^. All incoming and outgoing particles are taken to be on- mass-shell, 
kf = k2 = k^ = and P^ = Mj^^. In the static approximation each of the heavy quarks carries 
one half the mass and one half the four- momentum of the J/tp: 



Pc=Pc = =p; mc = Mji^/2 ; and 
From (Q) we find for the amplitude of Fig.2(a) 
' 1 
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Here g and e are the strong and electromagnetic couplings respectively, g = v^vra^, e = -v/Ittq, 
and Be is the magnitude of the charm quark charge in units of e. Charge conjugation invariance 
implies that the Feynman graphs are symmetric under reversion of the fermion flow. All six 
amplitudes contributing to the leading-order cross section are proportional to the same colour 
factor I/V3 6ij {T"-T^)ij = l/(2\/3) 5ab, where a and h are the colour indices of the gluonic quanta. 
Gauge invariance ensures that we may sum over the spins of the initial photon and over the gluon 
spins by employing the substitutions ^e^Ei, = —g^u- To sum over the J/^p spin states, we use 
X^Ej/^Ej/^ = —g'"^ + P^P'^ /Mj^^. Averaging over the initial photon/gluon helicities and colour 
the result for the cross section of the subprocess dH) may be written as 

da(0) _ I287r2 aa|e| . Iv^WP sj + tHj + uj 

Although the cross section ( |To|) is infrared finite it is not clear a priori in which kinematical 
region perturbative QCD can reliably be applied to J/^p photoproduction. Indeed, at large z the 
J/ip scatters more and more elastically and multiple soft gluon emission has to be considered. 
Similar effects might become important in the region where the transverse momentum p± of the 
J/ip tends towards zero. As discussed before, it is mandatory to require z <, 0.9 in order to eliminate 
contributions from elastic/diffractive production mechanisms. In Rcf . [p^] it has been argued that 
an additional cut on the transverse momentum of the J/ip, p± ^ 1 GeV, has to be applied to define 
the truly inelastic region. This conclusion was, however, based on a comparison of leading-order 
cross sections with experimental data that have later been found to be contaminated from coherent 



production in the small p± domain [22|. Accordingly, the region of applicability of perturbative 
QCD had not been clarified completely so far. As will become clear in Sec.^, the analysis of the 
next-to-leading order corrections restricts the kinematical domain where fixed-order calculations 
give a reliable description of inelastic J/ip photoproduction. We will observe that the perturbative 
QCD calculation is not under proper control in the singular boundary region z ^ 1 and p± — > 0, 
thereby indicating where multiple soft gluon emission becomes important. 



3 The NLO cross section 

Including higher-order QCD corrections to the partonic reaction (^) is expected to be essential 
for the theoretical description of inelastic J ftp photoproduction. Next-to-leading order corrections 
to open heavy flavour photoproduction have been calculated over the recent years mJH. They 



have been found to increase the normalization of the cross section significantly without strongly 
affecting the shape of the inclusive differential distributions. These results can however not directly 
be transferred to the case of bound state production as will become clear in Sec.^. It is thus 
important to investigate how the features of the lowest-order J/ijj photoproduction cross section 
are modified by radiative corrections and by including new production mechanisms which contribute 
in next-to-leading order. 

In this work we present a complete calculation of the higher-order perturbative QCD corrections 
to inelastic J/ip photoproduction. Results for total cross sections and the J/^p energy and transverse 
momentum spectrum will be discussed. The photon-parton reactions which contribute to the 



inclusive cross sections up to order aaf are 



7 + g {QQ)+g 0{aal), 0{aal) 

1 + 9 {QQ)+g + 9 0{aal 

1 + 9 {QQ)+q + q 0{aal 

i + q{q) ^ iQQ)+9 + qiq) 0{aal 



including virtual corrections to the leading-order process. We choose a renormalization and factor- 
ization scheme in which the massive particles are decoupled smoothly for momenta smaller than 
the heavy quark mass [^]. This implies that the heavy quark does not contribute to the evolution 
of the QCD coupling and of the structure functions. Furthermore, there are no contributing sub- 
processes initiated by an intrinsic heavy flavour coming directly from the structure function of the 
proton. All effects of the heavy quark are contained in the parton cross section. 

The calculation of the the next-to-leading order corrections will be outlined in this section. 
More details can be found in the Appendices. 



3.1 Virtual corrections 



The evaluation of the 0{as) corrections to inelastic J/tl^ photoproduction involves the calculation 
of the virtual cross section obtained from the interference term between the virtual and the Born 
amplitude. For the unrenormalized virtual cross section one finds in n = 4 — 2e dimensions 
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dtidui 



1 



1 



7r(47r) 



-2+e 



tiUi — Am?s 



-l)4(l-e)2 r(l-e) \^ 
x5(si +ti+ui+ 8m2) ^ 2Re(7W^7W^*) 



(11) 



where 1/{N'^ — 1) is the colour average factor for the gluon in the initial state and N denotes 
the number of colours. Photons and gluons have n — 2 spin degrees of freedom resulting in the 
spin average factor l/(n — 2)^ = 1/4(1 — e)^. The parameter /i has the dimensions of a mass and 
is introduced in order to compensate for the dimensionality of the gauge coupling constants in n 
dimensions. 

The Feynman gauge has been adopted to evaluate the 105 diagrams which contribute to the 
virtual amplitude. Charge conjugation invariance implies that the graphs are symmetric under 
reversion of the fermion flow. The ultraviolet (UV), infrared (IR), and the collinear or mass (M) 
singularities have been treated by using n dimensional regularization and show up as single and 
double pole terms of the type e~*(i = 1,2). We have refrained from adopting the seemingly 
simpler scheme of dimensional reduction which however gives rise to complications in the massive 



quark case as pointed out in Ref.|35|. The Feynman integrals containing loop momenta in the 
numerator have been reduced to a set of scalar integrals using an adapted version of the reduction 



program outlined in Ref.[36|. This program has been extended to treat n dimensional tensor 
integrals with linear dependent propagators in order to account for the IR- and M-singularities 
and the special kinematical situation which is encountered in the nonrelativistic approximation to 
J/ip photoproduction. The scalar integrals have been calculated by the Feynman parametrization 
technique and analytical results are listed in Appendix 

The UV-divergences which are contained in the fermion self energy graphs and the vertex 
corrections shown in Figj2|, are removed by renormalization of the heavy quark mass and the 



7 



'WWWWXA 




wvwwwv 



nnreireireinr*- 




(a) 



(b) 



'WWWWXA 



'WWVXAAAA 





(c) 



(d) 



Figure 2: Generic Feynman diagrams contributing to the virtual amplitude (part 1): self-energy 
corrections (a,b) [36 diagrams], abelian vertex corrections {jcc- and t^cc- vertices) (c) [18 diagrams], 
and non-abelian vertex corrections (d) [12 diagrams]. 



QCD coupling. The masses of the light quarks appearing in the fermion loops in Fig.^a have 
been neglected while the mass parameter of the heavy quark has been defined in the on-mass- 
shell scheme. The renormalization of the QCD coupling has been carried out in the extended MS 



scheme introduced in Ref.[34] and adopted in previous calculations of open heavy flavour production 
[ 37 , 3^,^] . In this scheme the effects of heavy flavours are smoothly decoupled for momenta much 
smaller than the heavy quark mass. The bare coupling has to be replaced according to 



ffba 




'jE + ln47r 




2 ( m 
/?0 + 3 In 



(12) 



with I3q = {UN — 2nf)/3. The total number of flavours including the heavy quark is denoted by 
nf and n\f = rif — 1 is the number of light quarks. The renormalization constant Zg is defined 
such that the contribution of the heavy-fermion loop in the gluon self energy graphs is cancelled 
in the renormalized cross section for small momenta flowing into the heavy-fermion loop. The 
UV-divergences arising from gluon or light fermion loops are removed according to the standard 



MS subtraction scheme [38]. For the renormalization scale evolution of the strong coupling one 
obtains from (12) 



5 In 



-"s(/^ij)/3o(nif) 



(13) 



The extended MS scheme thus implies that the heavy quark does not contribute to the evolution 
of the QCD coupling so that has to be evaluated using nif active flavours. 

Besides the self-energy diagrams and vertex corrections, 39 box graphs contribute to the virtual 
amplitude, which can be grouped in eight classes as shown in Fig.^. All diagrams falling into one 
particular class are related by exchange of photon or gluon momenta, adjustment of the colour 
factor and reversion of the fermion flow. These relations have been checked explicitly. 
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Figure 3: Generic Feynman diagrams contributing to the virtual amplitude (part 2): abelian 
four-point box graphs (a) [12 diagrams], abelian five-point box graphs (b) [6 diagrams], non-abelian 
four-point box graphs (c) [8 diagrams] and (d) [4 diagrams], non-abelian five-point box graphs (e) 
[4 diagrams] and (f) [2 diagrams], and diagrams involving four-gluon couplings (g) [2 diagrams] and 
(h) [1 diagram]. 
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The exchange of longitudinal gluons between the massive quarks in diagram Fig.^b leads to the 
Coulombic singularity ~ vr'^/u which can be isolated by introducing a small relative quark velocity 
V (see Appendix ^ for details). For the Coulomb-singular part of the virtual cross section we find 

M0r^'^''>(l + -CF- + ^C + O{a!)] , (14) 



where the colour factor is given hy Cp = {N'^ — 1)/{2N). The final state interaction in the colour- 
singlet channel is attractive and has to be interpreted as the Sommerfeld rescattering correction [^] 
which can be associated with the inter-quark potential of the bound state. The Coulomb-singular 
part of the virtual cross section is universal appearing in the next-to-leading order corrections to 
all production and decay processes involving 5- wave quarkonia. Following the standard path [pc| ], 
the corresponding contribution has to be factored out and mapped into the cc wave function: 



a 



\ vr V I 



l + -C + 0(a2) 
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(15) 



Only the exchange of transversal gluons contributes to the next-to-leading order expressions for the 
hard parton cross section. 



3.2 Real corrections 



The evaluation of the 0{aag) cross section requires the calculation of the gluon bremsstrahlung 
reaction 

7(^1) + 9{k2) ^ J/V'(2p) + giks) + gih) (16) 
and processes where the final-state gluon splits into light quark-antiquark pairs 



7(^1) +5(^:2) ^ J/^Pi2p)+q{k3)+qik4 



(17) 



The 48 Feynman diagrams which contribute to the amplitude can be obtained from the generic ones 
shown in Fig.^ by permutation of the photon and gluon lines. Spin and colour projection imply 
that the diagrams are invariant under reversion of the fermion flow. For the cross section of the 
two-gluon final states (|l6|), averaged over initial spins and colours, one obtains (see Appendix |B.1| ) 
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S3 = (A;3 + k^Y = si + ti + ui+ 8m^ and dl7„ = d6'i sin""'^ di (102 sin""'' 6*2. The angles Oi and 02 
which describe the orientation of the outgoing light partons are defined in Appendix |B.1| . A factor 
1/2! has to be included since there are two identical particles in the final state. For the cross section 
of the light-quark-antiquark-splitting reaction (|^) this factor has to be dropped. 
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Figure 4: Generic Feynman diagrams contributing to the gluon bremsstrahlung process (a-c) and 
the quark-antiquark sphtting reaction (d). 

The real-gluon cross section contains IR and M singularities so that the square of the amphtude 
had to be calculated in n dimensions up to order e^. For the sum of the gluon polarization we 
have used J2 ^f^^u = ~g^J.u and the unphysical longitudinal gluon polarizations have been removed 
by adding ghost contributions. 

The computation of the real cross section has been performed by adopting the phase space 
slicing method as outlined in Ref. ||4l]] and used in previous calculations of open heavy flavour 
production |35,33,42|. We have split the cross section into an infrared-collinear part (53 < A), 
which contains all singularities due to soft gluon emission and splitting of the final state gluon into 
gluon and light quark-antiquark pairs, and a hard-gluon part (53 > A). The cut-off parameter A 
is chosen such that it can be neglected with respect to mass terms like and the kinematical 
invariants si, ti, and ui. In the final answer the limit A ^ is carried out. The hard-gluon part 
(53 > A) contains single pole terms which are associated with initial state gluon radiation only. 
These collinear divergences have to be absorbed into the renormalization of the parton densities 



as outlined in Sec. 3^. To perform the integration over the orientation of the final state gluons or 
light quarks the matrix element squared has been decomposed into sums of terms which have at 
most two factors containing the dependence on the polar angle Oi and the azimuthal angle ^2- This 
decomposition is described in Appendix [B.l| . 

The infrared-collinear cross section (53 < A) for the two-gluon final states ( [T^ ) is obtained from 
the expression 
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For the light-quark-antiquark final state the symmetry factor 1/2! has to be dropped. The calcu- 



lation of the infrared-collinear cross section (20) is described in detail in Appendix B.2. Adding 



n 



the resulting expression ( [101 ,102) to the virtual correction leads to a cancellation of the infrared 
singularities. 



3.3 The photon-quark subprocess 

The cross section in next-to-leading order involves a new production mechanism where the photon 
is scattered off a light (anti-)quark from the proton 

lih) + q{q){k2) ^ J/^P{2p) + 5(^:3) + q{q){h) . (21) 

Because of spin and colour projection the J/V' particle can only be produced in the Bethe-Heitler 
reaction shown in Fig.^. The amplitude of the 'yq{q) subprocess does not depend on the electric 
charge of the light quarks. Averaging over spin and colour of the initial state particles the jq{q) 




Figure 5: Feynman diagrams contributing to the photon-(anti-)quark 

subprocess 79(9) — > J/il^gqiji). Additional graphs are obtained by reversing the arrows on the 
heavy quark lines. 



cross section is given by 

* dtidui ""^^ NA{l-e) 2r(l-2e) I iJ?s 



s^' J dflnY^lM^l^ . (22) 



The cross section ( p2[ ) contains collinear divergences ~ 1/e due to gluon emission from the initial- 
state light (anti-)quark which have to be removed by mass factorization. 

3.4 Mass factorization 

The collinear singularities contained in the hard-gluon bremsstrahlung reaction and the "yq{q) sub- 
process are universal and can be absorbed, as usual, into the renormalization of the parton densities. 
According to the factorization theorem |l| the parton cross section can be written as 



dtidu 



where s = xs, ui = xui and i = g,q{q). The reduced cross sections ddyg and d(T^g(g) defined by 
the above equation are free of collinear singularities and depend on the mass factorization scale 
Q^. This scale separates long and short distance effects and is a priori only determined to be 
of the order of the heavy quark mass m. The collinear singularities are contained in the splitting 



12 



functions Tgi of the incoming partons (gluons or light quarks). The sphtting functions depend on the 
mass factorization scale and further on the parameter /i^ which is an artefact of n-dimensional 
regularization. Up to leading order in they are given by p3| 



Tij{x, Q^, /u^, e) = 6ij5{l - x) + 



27T 



\Pi,{x) + fi,{x,Q^H^' 



(24) 



where the universal Altarelli-Parisi splitting kernels are denoted by Pij{x). The finite functions 
fij are completely arbitrary, different choices corresponding to different factorization schemes. Here 
we have adopted the MS-scheme corresponding to 



fff^ix, Q^,fx') = P,j{x) Le - ln47r + In ^ 



(25) 



For the reduced cross section of the 7-gluon process one finds from (23) 
d^o-^g (s,ti,Mi,(5^) _ d^cr^g (s,ti,Mi,//^,e) 



dtidui 



27r 







dx X Pgg{x) 



dtidul 

1 

1- 7-B - ln47r + In — ^ 

e /x^ 



d?a^g {xs, ti,xui) 



(26) 



The corresponding expression for the 7 — q{q) scattering reaction reads 
d^a^V(s,ii,iii,Q2) _ d^(T^y(s,ti,'Ui,/i2^g) 



dtidni 



as 



27r Jo 



dx X Pgq (x) 



dtidtii 

1 Q2 
1- 7-B - In 47r + In ^ 



d^a'^g{xs, ti, xui) 
dtidui 



(27) 



The Altarelli-Parisi kernel Pgg{x) of the gluon-splitting function has the form 



N 



Q{l-x-5)-2{ ^—^ + + x{l - x) 



1 — X 



+ 5(l-x) (21n5+y 



- nif 5(1 - x) 



(28) 



where the number of light flavours is denoted by nif. We have adopted the convention introduced 
in Ref. [pl| to regulate the pole at x = 1. The parameter 6 allows one to distinguish between soft 
(x > 1 — (5) and hard {x < 1 — 6) gluons and is related to the cut-off parameter A by 5 = A/(s-|-ui). 
Finally, the Altarelli-Parisi kernel Pgq appearing in (^) is given by 



Pggix) = Pgq{x) = Cp 



l + (l-x)^ 



(29) 



After mass factorization and cancellation of the infrared singularities between the virtual cor- 
rections and the contribution of soft gluon emission we obtain a finite expression for the 0{aa^) 
inelastic J/ip photoproduction cross section. An analytical result has been derived for the dou- 
ble differential one-particle-inclusive cross section da /dtidui. The corresponding expressions are 
however too long to be presented here. 
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4 The parton cross section 



The perturbative expansion of the total photon-parton cross section can be expressed in terms of 
scaling functions, 



aalel |v3(0)|^ 



ai^{s,ml) = — 1^—3- 



cS;)(7?)+4™Jc;;)(r/)+c;;)(r?)ln^ 



(30) 



where i = g,q,q denote the parton targets and (3o{ni{) = {UN — 2nif)/3. The scaling functions 
depend on the energy variable r] = s/Am^ — 1. c^g is the lowest-order contribution which scales 
~ r/~^ ~ Am^/s asymptotically. The cross section is put into a form in which the renormalization 
scale iiR and the factorization scale Q can be varied independently. The scaling functions c^i{Tj) 
are shown in Figs.^(a) and (b) for the parton cross sections integrated over z < zi where we have 
chosen zi = 0.9 as discussed before. [Note that the definition of z is the same at the nucleon and 
parton level since the momentum fraction x of the partons cancels in the ratio z = p ■ / p ■ k^.] 
In Fig.0 the scaling function of the gluon initiated parton process has been decomposed into a 
"virtual + soft" (V+S) piece and a "hard" (H) gluon-radiation piece. The In* A singularities of the 
(V+S) cross sections are mapped into (H), cancelling the equivalent logarithms in this contribution 
so that the limit A — > can safely be carried out. The nomenclature "hard" and "virtual + soft" 
is therefore a matter of definition, and negative values of c^^^ may occur in some regions of the 
parameter space. [In the range 0.2 ^ r] ^ 2 the hard gluon-radiation piece c^^'^^ as well as c^g 
differ from the curves in Ref.[^] by a few percent since the experimental cut z < 0.9 was not 
implemented properly in one term of Ref. ||2^ .] 



The following comments can be inferred from the figures, (i) The form of the hard-gluon 
radiation piece c^^-* resembles the corresponding scaling function in open-charm photoproduc- 



tion [33]. The logarithmic enhancement near threshold can be attributed to initial state gluon 
bremsstrahlung. The "virtual + soft" contribution for J/il^ production is, however, significantly 
more negative than for open-charm production. The destructive interference with the lowest-order 
amplitude is not unplausible though, as the momentum transfer of virtual gluons has a larger chance 
[in a quasi-classical approach] to scatter quarks out of the small phase-space element centered at 
Pc + Pc = Pj/^p than to scatter them from outside into this small element, (ii) While c^^^ and 
Cg^^^^^ scale asymptotically ~ 1/s, the hard coefficients Cg^^^ and Cq^ [as well as 0^57] approach 
plateaus for high energies, built-up by the fiavour excitation mechanism, (iii) The cross sections 
on the quark targets are more than one order of magnitude smaller than those on the gluon target, 
(iv) A more detailed presentation of the spectra would reveal that the perturbative analysis is not 
under proper control in the limit z ^ 1, as anticipated for this singular boundary region (see the 
detailed discussion in Sec.^). Outside the diffractive region, i.e. in the truly inelastic domain, the 
perturbation theory is well-behaved however. 



5 The photon-proton cross section 

The results for inelastic J/'0 production in photon-proton collisions 

7 + P ^ J/V' + X (31) 
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Figure 6: (a) Coefficients of the QCD corrected total inelastic \z < 0.9] cross section 7+5 J/ip+X 

in the physically relevant range of the scaling variable rj = s^p/Am? — 1; and (b) for 
1 + J/V' + X. 
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Figure 7: Coefficients of tlie QCD corrected total inelastic [z < 0.9] cross section 7 + 5 —> J/ip + X 
split into a hard piece and a virtual plus soft piece. 

are obtained from the partonic cross sections by convolution with the gluon and light-quark distri- 
butions ff in the proton 

da^P= ^ Jdxff{x,Q'')da^^ . (32) 

In the following we present a comprehensive analysis of total cross sections and differential distri- 
butions for the energy range of the fixed-target experiments and for inelastic J/tp photoproduction 
at HERA. 



5.1 The energy range of the fixed target experiments 

Inelastic J/ip photoproduction has been measured in fixed-target experiments 12^,21] at photon 
energies near E-y = 100 GeV, corresponding to invariant energies of about ^/s■yp 14 GeV. Before 
comparing the theoretical predictions with the experimental data we will examine the effect of the 
next-to-leading order corrections in some detail. 

In Figs.|8|and^the J/ip energy spectrum da/dz and the J/x/j transverse momentum distribution 
da/dpj_ are shown at an initial photon energy of E^y = 100 GeV. The GRV parametrizations of the 



parton densities |Q have been adopted. They are particularly suited to characterize the magnitude 
of the radiative corrections since they allow one to compare the results for the Born cross section 
folded with leading-order parton densities, with the cross sections consistently evaluated for parton 
cross sections and parton densities in next-to-leading order. As the average momentum fraction 
of the partons <a;>~ 0.1 is moderate, the curves are not sensitive to the parametrization in the 
small-x region. The renormalization scale has been identified with the factorization scale and set to 

(5) 

Hr = Q = Mj/^. For as the two-loop formula is used with nif active flavours and = 215 MeV, 
corresponding to the average fit value in Ref . . Since the cross section depends strongly on the 
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Figure 9: Transverse momentum distribution da/dp]_ at an initial photon energy of = 100 GeV 
integrated in the inelastic region z < 0.9. 
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QCD coupling, we adopt this measured value, thus allowing for a slight inconsistency to the extend 
that the GRV fits are based on a marginally lower value of Og (cf. Ref.||4^). In next-to- leading 
order, the wave-function at the origin is related to the leptonic J width by 

16 Os 
TV 



1 



(33) 



with only transverse gluon corrections taken into account explicitly 
3.097 GeV |6| and rric = Mj/^/2. 



We use Tp 



5.26 keV, 



M 



Prom Fig.g one can infer that the perturbative QCD analysis is not under proper control in 
the limit z ^ 1, as anticipated for this singular boundary region. If we restrict the analysis to 
the inelastic domain z <, 0.9 the perturbative expansion is well-behaved however and the next- 
to- leading order corrections do not strongly affect the shape of the distributions. The i^-factor, 
K = (Tnlo/clO) is nearly independent of z and p± in the inelastic region z <^ 0.9. Its magnitude 
turns out to be ~ 2.0 with one part ~ 1.73 due to the QCD radiative corrections of the leptonic 
J/ip width and a second part ~ 1.2 due to the dynamical QCD corrections. The slope of the 
transverse momentum distribution, da/dp^ oc exp{—bp'j_), is predicted to be 6 0.6 GeV~^, in 
good agreement with the experimental value b = 0.62 it 0.2 GeV~^ [20|. 

The scale dependence of the theoretical prediction is reduced considerably when higher-order 
QCD corrections are included. This is demonstrated in Fig.|l^ where we compare the scale depen- 
dence of the leading and next-to-leading order total cross section in the inelastic region z < 0.9. 
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Figure 10: Dependence of the total cross section j + P J/ip + X on the renormaliza- 
tion/factorization scale Q at an initial photon energy of = 100 GeV. 



For the sake of simplicity, the renormalization scale has been identified with the factorization scale, 
/^i?, = Q- While the ratio of the cross sections in leading order for = : (2m^) : ^j/^ is given 
by 1.7 : 1.3 : 1, it is much closer to unity, 0.7 : 1.1 : 1, in the next-to-leading order calculation. The 
cross section runs through a maximum |49| near ~ 2m^ with broad width, the origin of the 
stable behaviour in Q. 
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The next-to-leading order predictions can be confronted with photoproduction data of the fixed 
target experiments. For a meaningful comparison, the theoretical uncertainties due to variation 
of the charm quark mass and the strong coupling have to be taken into account properly. In 
the static approximation the choice rric = is required for a consistent description of the 

bound state formation. A smaller mass value, however, is favoured to describe the charm quark 
creation in the hard scattering process [^]. The value of the heavy quark mass in the short 
distance amplitude is the main parameter controlling the normalization of the cross section. It 
is therefore appropriate to adopt charm masses below rric = thus allowing for a slight 

correction in the bound state formation. Leading-order analyses of J/tp production that go beyond 
the static limit and incorporate a non-vanishing binding energy, find an effective charm mass value 
of rric = 1.43 GeV [p^ j, in fairly good agreement with potential model calculations [^]. In order to 
demonstrate the uncertainty due to the variation of the charm quark mass, the strong coupling and 
the renormalization/factorization scale the results will be shown for (i) rric = -^j/i/)/2 ~ 1.55 GeV 
215 MeV and = fi'j^ = Mj, , (as used in the previous figures) and (ii) rric = 1-4 GeV 



,(5) 



with A.^g 

With Ag 



= 300 MeV (corresponding to the 1 a upper boundary of the error band in Ref.|4|]) and 
2mt 



In Fig. 11 we confront the leading and next-to- leading order calculations with the J/ip energy 
spectra measured at photon energies near = 100 GeV HJH. It is clear from Fig.|l| that 
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Figure 11: Energy spectrum da/dz, at the initial photon energy = 100 GeV compared with the 



photoproduction data |20| , |21[| . 



the variation of the charm mass and the strong coupling does not strongly affect the shape of 
the distribution but only results in some uncertainty concerning the overall normalization. In a 
systematic expansion one may finally add the relativistic corrections as estimated in Ref . pSj . 



The dependence of the total cross section 7+P J/ip+X on the photon energy E^ is presented 



in Fig. 12, again for the two choices of parameters (i) and (ii) as defined above, together with the 
photoproduction data ||20|j2l[| . From the curves shown in Fig.|l2| we deduce that the QCD corrections 
are large at moderate photon energies, but decrease with increasing energies, a consequence of the 
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negative dip in the c^^^ scaling function of FigJ 
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Figure 12: The total cross section [z < 0.9] as a function of the initial photon energy. 



Two conclusions can be drawn from the comparison of the next-to-leading order results with 
the experimental data, (i) The J/ip energy dependence da/dz and the slope of the transverse 
momentum distribution da/dpj_ are adequately accounted for by the theoretical prediction in the 
inelastic region, (ii) The absolute normalization of the cross section is somewhat less certain. 
However, taking into account the theoretical uncertainty due to variation of the charm quark mass 
and the strong coupling and allowing for higher-twist uncertainties of order (A/ttIc)'^ ^ 30% for 
A; > 1, we conclude that the normalization too appears to be under semi-quantitative control. 



5.2 Inelastic J/ip photoproduction at HERA 

The production of J/V' particles in high energy ep collisions at HERA is dominated by photopro- 
duction events where the electron is scattered by a small angle producing photons of almost zero 
virtuality. The measurements at HERA provide information on the dynamics of inelastic J/ ip pho- 
toproduction in a kinematical region very different from that available at fixed target experiments. 
The 7p centre of mass energies accessible at HERA are in the range 30 GeV ^ \/sjp ^ 200 GeV, cor- 
responding to initial photon energies in a fixed-target experiment of 450 GeV ^ E.y <, 20, 000 GeV. 

To begin with, we discuss the J/jp energy spectrum da/dz and the J ftp transverse momentum 
distribution dc/dp^ at a typical HERA energy of ^/s^p = 100 GeV. The parameters have been 
chosen as in the corresponding figures of Sec. 5.1. From Figs.^ and 14 one can conclude that the 



next-to-leading order corrections are dominated by strong negative contributions in the limit z ^ 1 
and p± 0. This behaviour which has already been observed in the low energy region of the 
fixed-target experiments is much more pronounced in the high energy range at HERA. Even if 
the analysis is restricted to the region z < 0.8 we still find that the fixed-order perturbative QCD 



calculation is not under proper control for p± 0. This can be inferred from Fig.|l^ where the 
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transverse momentum spectrum is shown integrated in the range z < 0.9 and z < 0.8, respectively. 
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Figure 13: Energy distribution dcr/dz at the photon-proton centre of mass energy ^Js^y =100 GeV 
integrated in the full p_L range and in the restricted range p± > 1 GeV. 
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Figure 14: Transverse momentum distribution dcr/dp^ at the photon-proton centre of mass energy 
^ s^p = 100 GeV integrated in the region z < 0.9 and z < 0.8. 

For p_L ^ 0.5 GeV the results of this calculation obviously require missing contributions from even 
higher orders in the perturbative expansion. No reliable prediction can be made in the small p± 
and large z domain without resummation of large logarithmic corrections caused by multiple gluon 
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emission. It is therefore appropriate to exclude the region z ^ 1 and p_l ^ from the analysis. 
In the following we will present the results in two kinematic domains, (I) z < 0.9, which is the 
minimal restriction in order to eliminate elastic/diffractive contributions, and (II) z < 0.8 and 
p± > 1 GeV, which was found to be the region where fixed-order perturbation theory allows for 
a reliable prediction in the HERA energy range. The next-to- leading order results for the J/V' 
energy spectrum are shown in Fig.^ integrated in the full p± range and in the restricted range 
p± > 1 GeV. No singular behaviour is observed for the latter curve even in the limit z ^ 1. 



In Fig, 15 we compare the scale dependence of the leading order and next-to-leading order 

100 GeV. The renormalization scale has been 
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Figure 15: Dependence of the total cross section j + P J/ip + X on the renormaliza- 



tion/factorization scale Q at the photon-proton centre of mass energy y/s 



IP 



100 GeV. The 



results are shown in two kinematic domains: (I) z < 0.9; (II) z < 0.8 and p± > 1 GeV. 



identified with the factorization scale, = Q. The results are shown in the two kinematic domains 
(I) z < 0.9 and (II) z < 0.8 with p± > 1 GeV, as discussed above. From Figjl^ one can infer that 
the next-to-leading order result is insensitive to scale variations in an appreciable range near the 
reference scale /i^ = = ^j/^p- For scales below ii\ = ^ no stable prediction 

is possible. In contrast to the low energy region, the cross section in the HERA energy range 
does not exhibit a point of minimal scale sensitivity. In the BLM scheme |5^ we find a value of 
~ This scale is significantly larger than the corresponding BLM value for J/^/j decays. 

The typical kinematical energy scale is not set any more by the small gluon energy in the J/'i/' 
decay but rather by the typical initial-state parton energies. 

We will now present our final predictions for differential distributions and total cross sections for 
inelastic J/ij: photoproduction at HERA. In a systematic expansion we have added the relativistic 
corrections [^] which enhance the large z and small region and thereby increase the total cross 
section integrated in the range z < 0.9 by ~ 10%. The inclusion of relativistic corrections does 
not change the results obtained in the more restricted domain z < 0.8 and p± > 1 GeV. In order 
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to demonstrate the theoretical uncertainty the results are shown for (i) = ~ 1.55 GeV 

215 MeV and = fij^ = Mj^^ (as used in the previous figures) and (ii) rUc = 1.4 GeV 
300 MeV and = 



with A§ 
with Ag 



2 m^, as discussed in Sec. |5.l| . In Figjl^ we plot the transverse 



momentum distribution at the photon-proton centre of mass energy ^/s 



IP 



100 GeV integrated 



in the region z < 0.9 and z < 0.8. As could already be inferred from Fig, 14, the inclusion of the 
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Figure 16: Transverse momentum distribution da/dp'^ at the photon-proton centre of mass energy 
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100 GeV integrated in the region z < 0.9 and z < 0.8. 



next-to-leading order corrections increases the cross section in the range p_l ^ 1 GeV and results 
in a hardening of the distribution. For the slope, da/dp'^ oc exp(— we predict b ~ 0.3 GeV~^. 
The J/ip energy distribution is shown in Fig.^ integrated in the full p± range and in the restricted 
range p± > 1 GeV, again for two choices of the charm quark mass and the strong coupling. In 
Figs.^(a) and (b) we plot the total cross section as a function of the photon-proton centre of mass 
energy in the HERA range. The results are shown in the two kinematic domains (I) z < 0.9 and 
(II) z < 0.8 with p_L > 1 GeV. For domain (I), Fig.|l8|(a), the /C-factor is ~ 0.75, a consequence of 
the strong negative contribution present in the region p_l — > 0. In the more restricted kinematic 



domain (II), Fig. 18(b), the next-to-leading order corrections significantly increase the cross section. 

,(5) 



For rrir = IA GeV with A^ 



1.7. 



300 MeV and Q"^ = p% = 2ml, we find K = ctnlo/o-lo 

An estimate of the cross section for inelastic photoproduction of ip' particles can be obtained 
from the results presented here by replacing the leptonic decay width and multiplying with a phase 
space correction factor 



a(7P ^ ^' X) 



l/Ax a{-fP ^ J/ip X) . 



(34) 



The estimate (34) should be considered as a lower bound since it is based on a purely static 
approach. In the derivation of the phase space suppression factor (Mj/^/M^/)^ it is assumed that 
the effective charm masses in the short distance amplitudes scale like the corresponding ip' and J/V' 
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Figure 17: Energy distribution dcr/dz at the photon-proton centre of mass energy ^s^p = 100 GeV 
integrated in the fuh range and in the restricted range p_l > 1 GeV. 



masses. This difference might be reduced by including relativistic corrections which are expected 
to be significantly larger for the "0' state than for the J/'i/' |51]. A measurement of the inelastic '0' 
photoproduction cross section at HERA would clearly help to understand the impact of relativistic 
corrections on charmonium production. 

The production of T bottomonium bound states is suppressed, compared with J/^/j states, by 
a factor of about 300 at HERA, a consequence of the smaller bottom electric charge and the phase 
space reduction by the large h mass. 

Since the momentum fraction of the partons at HERA energies is small, the cross sections 
presented above are sensitive to the parametrization of the gluon distribution in the small-x region 
<x>~ 0.003. In Figs.|T9|(a) and (b) we compare the next-to-leading order predictions for different 
parametrizations of the gluon distribution in the proton with first results measured at HERA 
HH]. We have used = 1.4 GeV with A^ = 300 MeV and Q'^ = ii\ = 2ml, favoured 
by the low energy data. The data samples [ p3| , |54[ contain background events from production of 
V'' states with subsequent J decay |18|. This contribution has been included in the theoretical 
prediction, Figs.[l9|(a) and (b), in order to allow for a meaningful comparison. According to (]3^ ) and 
the measured branching ratio BR(^' J/iIjX) = 57 % |Q the cascade production from •0' states 
is conservatively expected to increase the cross section by ~ 15%. For the parametrizations of the 
parton densities in the proton we have chosen four sets that are compatible with the recent HERA 
measurements of the proton structure functions |^5|: the GRV parametrization ||45[| , which has been 
adopted in all previous figures, the sets MRS(G), MRS(A') |Q and the CTEQ3 parametrization 
[^7[ . The difference in the small x behaviour of the gluon densities, xg ~ x~^ where Agrv ~ 0.3—0.4, 
Amrs(g) ~ 0-4, Amrs(A') ~ 0-2 and ActeQ3M ~ 0.3, results in different normalizations and, to a 
smaller extend, in different shapes of the cross section as a function of the photon-proton centre of 
mass energy. The shape of the differential distributions in z and p±_ is not very sensitive to changes 
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Figure 18: The total cross section for inelastic J/V' photoproduction 7 + P J/tp+X as a function 
of the photon-proton centre of mass energy in the HERA energy range. The results are shown in 
two kinematic domains: (a) z < 0.9; (b) z < 0.8 and p± > 1 GeV. 
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Figure 19: The total cross section as a function of the photon-proton centre of mass energy for 
different parametrizations of the gluon distribution of the proton compared with prehminary data 
from HI |53| and ZEUS [54|. The results are shown in two kinematic domains: (a) z < 0.9; (b) 
z < 0.8 and p± > 1 GeV. 
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in the proton structure function. The theoretical result for the total cross section integrated in the 
kinematic domain z < 0.9 on average underestimates the experimentally observed production rate, 
Fig.pj|(a). This could be anticipated since the perturbative analysis is not under proper control 
in the limit p± ^ and the next-to-leading order corrections are dominated by strong negative 
contributions for p± <, 0.5 GeV, as shown in Fig. 14. Following the arguments above, it is more 
adequate to compare theory and experiment in the kinematic domain z < 0.8 and p± > 1 GeV, 
where fixed-order perturbation theory allows for a reliable prediction in the HERA energy range. 
This comparison is shown in Fig.|l9|(b). It can be inferred from the plot that the next-to- leading 
order result not only accounts for the energy dependence of the cross section but also for the 
overall normalization. The sensitivity of the prediction to the gluon distribution in the proton is, 
however, not very distinctive in the more restricted domain z < 0.8 and p± > 1 GeV. In particular 
the MRS(G) and GRV parton densities lead to almost identical results over the whole kinematical 
range accessible at HERA. A detailed analysis reveals that the size of the QCD corrections increases 
when adopting parton densities with flatter gluons. The sensitivity to different gluon distributions 
is thus reduced in next-to- leading order as compared to the leading-order result, in particular when 
choosing a small charm mass and a large value for the strong coupling. Parametrizations with 
extremely flat gluons like MRS(DO') |^8| are clearly disfavoured by the recent HERA measurements 



of the proton structure function |55| and do not allow for a reliable prediction in the high energy 
region. For the parameters adopted in Fig. 19(b), the MRS(DO') distribution leads to next-to- 
leading order results not very different from those obtained with the MRS(A') parametrization. The 
corresponding X-factors are however uncomfortably large, X ~ 4, casting doubts on the reliability 
of the perturbative expansion as obtained by using flat gluon distributions. If parton distributions 
with steep gluon densities are adopted, the next-to-leading order cross section is well-behaved and 
gives an adequate description of the experimental data, as demonstrated in Fig.pj|. 



6 Conclusion 



We have presented a complete calculation of the higher-order perturbative QCD corrections to 
inelastic photoproduction of J/V' particles. In the energy range of the fixed target experiments, 
Ery ~ 100 GeV, including the next-to-leading order corrections reduces the scale dependence of 
the theoretical prediction and increases the cross section by about a factor of two. A comparison 
with photoproduction data of fixed-target experiments reveals that the J/^p energy dependence and 
the slope of the transverse momentum distribution are adequately accounted for by the theoretical 
prediction in the inelastic region z ^ 0.9. Taking into account the uncertainty due to variation 
of the charm quark mass and the strong coupling and allowing for higher-twist effects of order 
(A/mc) <^ 30%, we conclude that the normalization too appears to be under semi-quantitative 
control. In the high energy range at HERA, a detailed analysis of the spectra has shown that 
the perturbative calculation is not well-behaved in the limit p± — > 0. No reliable prediction can 
be made in this singular boundary region without resummation of large logarithmic corrections 
caused by multiple gluon emission. First experimental results from HERA indeed indicate that 
the production rate, obtained in the full p± range, is on average underestimated by the theoretical 
predictions. If the small p± region is excluded from the analysis, the next-to-leading order result not 
only accounts for the energy dependence of the cross section but also for the overall normalization. 
The results seem to favour a gluon density in the proton rising toward low x, consistent with recent 
measurements of the proton structure functions. Higher-twist effects must be included to improve 
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the quality of the theoretical analysis further. 
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APPENDICES 



A The scalar integrals 



Analytical results for the scalar integrals which emerge from the Passarino-Veltman reduction of 
the virtual amplitude are listed in this Appendix. The exchange of longitudinal gluons between 
the massive quarks in diagram Fig.^a leads to a Coulomb singularity ~ Tr^/ii, which appears in the 
evaluation of the corresponding loop integrals. The singularity has been isolated by introducing 
a small quark velocity f , as discussed in detail below. First, we will give the expressions for the 
Coulomb-finite integrals. The results are not analytically continued into the physical region 



5,51+ 4m^ > 4m^, 



t, ti + W < 0, 



u , til + Am < 



(35) 



with si, ti and ui as defined in (0). Therefore, the expressions for all permutations of the photon 
and gluon momenta can be obtained from the formulae listed below by interchanging ti ui, 
si ^ ti, or si ui. Since any imaginary part of the integrals will disappear in the final result, 
only the real parts are given. The four-momenta are related by ki + k2 = 2p + and all particles 
are taken to be on-mass-shell, /cf = /c^ = A^l = and 



m 



The scalar one-point function is given by 

d"g 



A{m^) = n 



4— n 



where 



(2vr) 



1 



n q2 



iCf m^ 



+ 1 



167r2 



-e(7B-ln47r) ( 



m^ 



(36) 



(37) 



and n = 4 — 2e. 

— The scalar two-point function is defined by 



B{p,mi,m2) = n 



A-n 



(38) 



{2tt)^ [g2 — mf] [{q + pY — ml^ 
The following types of two-point functions appear in the calculation of the virtual amplitude: 

S(A;i,0,0) = (39) 



28 



5(^2-^3,0,0) = iC, 

B{p,0,m) = B{2p,m,m) = iC, 
B{p — ki,0,m) = iCe 
B{ki + k2,m,m) = iC, 



--In (—^1+2 



1 

- + 2 
e 



- + 2 — - In , 

e ti + 2m2 V 2m2 



1 



+ 2 + /31n 



V1 + /3; 



Here and below we have used the shorthand notation /3 = = ^1 — Am^ /s. 
— The scalar three-point function is defined by 

C{pi,p2, mi, 1712,1713) = 
A-n f d"(7 



-n f 

y (27r)" 



1 



(27r)" [g2 - ml] [{q +pif- ml] [{q + p2f - ml] 
Eight different types of three-point functions appear in the virtual amplitude: 



C(A;2,fe3, 0,0,0) = iC, 
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4m2 - s(l - 13) 
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C{ki, ki + k2,m, m, m) 
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1-/3 



where C(2) = 

— The scalar four-point function is defined by 

D{pi,p2,P3,mi,m2,m3,mi) = 
A-n f 



1 



(27r)" [g2 - J7i2j [(^g ^ p^)2 _ J(^g ^ _ j^2j _)_ ^^-^2 _ ^ 

There are three different types of four-point functions: 
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where A = A(si, ii) 



Sl ti + 4m2 



(55) 



(56) 



y Sl + 4m^ ti 

The five-point functions appearing in the calculation of the Feynman amplitudes shown in Fig.^(e,f ) 
can be rewritten as linear combinations of the four-point integrals listed above. 

In order to evaluate the Coulomb-singular diagrams shown in Fig.^(b) three five-point tensor 
integrals of the type 

Eo,a,af3,af3'y,af3'y5{Pc,Pc - ki,Pc- ki - k2, -Pc, 0, m, m, m, m) = 

4_„ f d"g I,qa,qaqi3,qaqi3q'y,qaqi3q-yq5 



(27r)" [(q + p,)2 _ m^] [(g +p^-k^)2_ ^2] 

1 



[{q +Pc- ki- k2Y - m2] [{q - p^)^ - rr?] 



(57) 



have to be calculated. The two additional integrals can be obtained from (57) by exchanging the 
photon and gluon momenta. The Coulomb singularity has been isolated by introducing a small 



quark velocity v and the momenta p^. and fulfill the relations Pc + Pc = and \pc — Pc\ = mv. 
Since all loop integrals which have at least one power of the loop momentum q in the numerator 
are Coulomb-finite, the tensor reduction of (|57|) can be performed at v = 0. The vr^/t; singularity 
is thus contained only in the scalar five-point integrals Eq and in four- and three-point functions 
of the type D{pc,Pc — ki, —pc,0,m,m,m) and C{pc-,—Pc-,^ifn,m) which emerge from the tensor 
reduction. The integration of the three-point function is straightforward and we obtain in the static 
limit V 

1 



C{pc,-Pc,0,m,m) 



-ia 
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1 ^ 

e V 



(58) 



Alternatively, the Coulomb singularity (as well as the infrared pole) can be isolated by introducing 
a small gluon mass A. In the limit A — > one finds 



C{pc, -Pc, 0, m, m) = -iC^ 



1 
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W (4) 


+ — m - 2 
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(59) 



which shows the correspondence between the different regularization schemes (cf. Ref . j59|| ) : 

(IR-singularity) 

(Coulomb-singularity) . (60) 

The singular piece of the four- and five point functions is determined by the infrared part of the 
integration region, g — > 0, and can thus be extracted by neglecting the loop momentum in the 
infrared finite propagators. In the limit \pc — Pc\ = "lu — > one finds for example 
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(61) 



and similar expressions for the other four- and five-point functions. No terms linear in v appear 
in the relativistic expansion of the coefficients multiplying the Coulomb-singular integrals. The 
corresponding expressions can thus be evaluated in the static limit v = from the start. 

For the finite part of the four-point integral we obtain 

(PcPc- ki,-pc,0,m,m,m) = 

D{Pc,Pc - ki , -pc, 0, m, m, m) 
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Similarly, the finite part of the five-point integral is given by 

E''^'^\pc,Pc - ki,pc- ki - k2,Pc,0,m,m,m,m) = 

E{PciPc - ki,Pc - ki- k2,Pc,0,m,m,m,m) 
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m^ 



Cipc,-Pc,0,m,m) 



(63) 



[{Pc - kiY - m'^][{pc -ki- k2y 

The difference = E — E^^^^ is infrared- and Coulomb-finite and a complete analytical result 

has been obtained. Since the corresponding expression is rather extensive we will only outline the 
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basic steps of the calculation. By using the Feynman parametrization technique one obtains 

i f fi^y r(3 + e) ^ 
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(4vr)n/2 



where 
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2ki ■ k2 



Substituting 

X = X + y — 1 , y = w , z = z — w , w = x — z , 
we arrive at an expression which does not depend on the integration variable x: 



dx 



2-2x rx 



dy dz dw {1 — y) — tyz — uyw — sw{l — z) 

10 Jl-x Jo Jo ^ J 

Interchanging the order of integration and evaluating two integrations, we find 
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The singularity is contained in the expression 

= 7^ t d2 /' dy {y - 1) f{e, s, t, u, z, y) 

l + € Jo Jo 



with 
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The integration of is straightforward and we obtain 
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The remaining difference 

= dz J\y[{y-z){fie,s,i,u,z,y)+2f{e,s,i,u,z,2y)} 

-{y- f'^{(^,s,i,u,z,y) 
is finite and can be evaluated in the limit e — > by using elementary integration techniques. 
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B The three-particle final states 



In this Appendix we outline the kinematics for the processes involving three-particle final states 
and give explicit results for the infrared-collinear cross section. 



B.l The phase space integration 



ti = {2p - kif - An? u' = {ki - ^2^^ 



For the calculation of the two-to-three-body processes (||), dl^) and (^) we closely follow Ref.[||] 
and introduce the following 10 kinematical invariants: 

s = {ki + k2f = si+ An? ta = (fcg - fcs)^ 

53 = ih + kif U6 = {k^-kif 

54 = {2p + kif - An? Ml = {2p - fca)^ - An? (75) 

55 = {2p + k^f-An? t' = {ki-kif 

where ki + k2 = 2p + k^ + ki^. The invariants s, ti and ui have already been defined in (0). Since 
we are dealing with a three-particle final state only five of the invariants are linearly independent. 
In the centre of mass frame of the outgoing light particles (gluons and/or light (anti-)quarks), the 
momenta are given by 

ki = (wi, . . . ,0,0,a;i) 
Pj/i> = '^P = {Ej/^,...,0,pj/^sm-4),pj/^cosi^) 

k^ = {ujs, ... ,—u)3sm9ism92,—uj3sin9icos92,—u;3COs9i) 

ki = (a;4, . . . , 0^3 sin 01 sin 6*2, sin^i cos 6*2, cos ^i) 

k2 = {u;2,...,0,pj/^sm^,pj/^cos^p-u;i) (76) 

The coordinate axes can be chosen in such a way that the n — A unspecified angular components 
do not contribute to the matrix element squared. From four-momentum conservation and the 
on-mass-shell constraints one can derive the identities 

•51 + ^1 + 4m^ si + til + An? ^/sz 
u^i = z-7= , W2 = —j= , u:^=Ui = —- 

^ Si - S3 r=r2 -r^ I ^1 -^2 

''■^/^= V^' W = ^/^./^-4m^ cos^P= 2^^^^^^ . (77) 
For the cross section one has 

^ = 2^.''-'''Ys^' ' J J2^ J J2^ J J2^ 

x5+(fci) 5+(fc|) (27r)" 5"(A:i + k2 - Pj,^ - k^ - k^) jA^^p , (78) 

Kc and Kg denoting the colour and spin averaging factors. The factor 1/2! has to be included for 
the two-gluon final states because of Bose symmetry. It is convenient to rewrite the cross section 
according to 

1 ,.2(4-n) f. 

a = KcKs- -^^2^ j d"Pj/^ - Aw?) 5^{ki + k2 - Pj/^ - 

1 
2! 



xi /d"fc3 /d"A:4<5+(fci)5+(A:|)5"(C-A:3-fc4) lA^^l' • (79) 



The integral over the phase space of the two outgoing light partons is evaluated in the corresponding 
centre of mass frame: 

^ J d^ks J d"A;4 6+{kl) 6+{kl) <5"(e - ks - k,) \Mr\^ = (80) 

^(n-4)/2 ^(„_4)/2 1 Ilizi) 1 r T |_^ij|2 _ 

4 r(n - 3) 2 7o Jo 

The remaining integrations are performed in the centre of mass frame of the initial state particles 
with momenta given by 



ki = ^/i/2(l,..., 0,0,1) 

k2 = V^/2(l,..., 0,0,-1) 

Pj/^ = {Ej/^,. . . ,0,pj/^smx,pj/^cosx) ■ (81) 
One finds for the cross section 

1 (4^)-"/u2(4-") f , 
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^ = ^^^^7 ^'rr^^.^ /di^^/^(i?^/^-4m^)-2- dxsin'^-^^ 

1 
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After changing the integration variables {Ej/^,x) {ti,ui) using 

we finally obtain 

2 dV ^ ^ (4^)->^-" Ai^i-4m^. y"-^^/^ 
^ dtidui ^ ^ 2r(n-3) 1^ /x2s j 

^^(n-4)/2l r^^^sinn-s^ /"^02sin"-^e2 lA^^I^ • (84) 
2 JO JO 

The kinematical limits on ti,ui can be deduced from the constraints —1 < cosx ^ 1 and S3 > 0: 

^THj^ s 

-s <ti< -4m^ - si-h- 8m^ < tii < . (85) 

To perform the angular integrations the matrix element squared has been decomposed into 
sums of terms which have at most two factors containing the dependence on the polar angle 61 
and the azimuthal angle ^2- This partial fractioning exploits the kinematical relations between the 
invariants defined in (|75|): 
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By using identities as e.g. 
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(86) 



(87) 



the angular integration can be reduced to the evaluation of standard integrals of the form 



j(k,i) 



r^Oi sin"-3 Oi rd92 sin^ 
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U.U2 air"" ^ ' 
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The results for these integrals can be found in Ref.[35 



B.2 The infrared-collinear cross section 

The calculation of the infrared-collinear cross section ( pOD will be outlined in this Appendix. The 
resulting expressions contain all divergences due to soft gluon emission and splitting of the final 
state gluon into gluon and light quark- antiquark pairs. 

We begin by considering the infrared region of phase space where one of the final state gluon 
momenta k^^^ becomes soft. In the limit A;4 — > the amplitude of the real gluon emission factorizes 
into an eikonal factor multiplying the Born amplitude: 



and analogous for A^s — > 
where = N{N'^ - 1). 



rr2(T«T'')=l 



^k^-k^ b 2 (90) 



^3^0 S3fg 2 



In addition to the soft gluon divergences, the amplitudes for the 2^3 processes contain terms 
which are singular when the two outgoing light partons become collinear, A;3||A:4, i.e. cos034 = 
COS0 = 1. For the calculation of the collinear matrix element we choose the four-momenta in the 
centre of mass frame of the incoming partons according to 

ki = — (l,...,sinx,0,cosx) 
v/i 

h = — (l,...,-sinx,0, -cosx) 
ks = a;3(l,..., 0,0,1) 
^4 = a;4 (1, ... ,sin0cos(/), sin sine/), cos^) 
P = 2p = (£',..., —a;4 sin cos 0, —a;4 sin 6* sine/), —1^4 cos 6* — a;3) . (91) 

Hence, in the limit 9^0 the condition S3 < A is given by 

S3 = 2 W3a;4 (1 - cos 6*) 6^36^4 6*^ < A . (92) 

When integrated over the solid angle 

an oc d0 sin"-^ 9 w e''~^d9 , (93) 

the I/S3 ~ 1/9'^ terms in the matrix element squared lead to logarithmic singularities in the limit 
0^0. The pieces with additional powers of 9 in the numerator only give contributions of the 
©(A) which vanish in the limit A ^ 0. Accordingly, all invariants in the expressions multiplying 



the I/S3 propagator can be replaced by their values at 6 = 0: 
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The other invariants remain unaltered. The splitting parameter z is defined as 

UJ4 



z{ti + 4m^) 
(l-z)(ti+4m2) 
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(94) 
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W3 + ^4 

and should not be confused with the J/ip energy variable defined in (^). The parts of the matrix 
element squared which are proportional to I/S3 ~ 1/6*^ have to be expanded up to the 0{6'^). Gauge 
invariance ensures that in the resulting expression all terms ~ 1 cancel. The contributions ~ 6 jQ'^ 
are proportional to cos(/> and vanish after the azimuthal integration. The parts ~ 9'^/ 6^ finally 
contribute to the collinear matrix element and lead to logarithmic singularities when integrated 
over the solid angle (|93|). 



For the process ^ J/tp + gg we obtain the following result for the matrix element squared 
in the limit 6* — > 0: 
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where B{x; y) = r(2;)r(y)/r(x + y) is the Beta-function. The gluon-gluon splitting function Pgg is 
given by 
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As mentioned before, the terms of the order 1/9 which have not been included in ( p6| ) only give 
contributions of the 0{A) and thus vanish in the limit A ^ 0. 



The infrared-collinear region of phase space (^) is illustrated in the ^;-0-plane in Fig.2C. For 
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Figure 20: The infrared-col- 
linear region of phase space 
in the z-9-plane. The in- 
frared parameter 5 and the 
collinear parameter e are ob- 
tained from S3 = A with 
S3 ~ 2; (1 — z) (1 — cos 6*) and 
vanish when A ^ 0. 



A;4 ^ (z ^ 0) and /c3 — > (z — > 1), the matrix element squared can be approximated by the 
expressions (|8^) and (|90|), respectively, while for ^ the collinear limit (p^) has to be used. It is 



clear from (|89|), (9C) and ( |96[ ) that in the infrared-cohinear region, i.e. in the upper left and right 
corner in Fig.2C, both approximations are identical. 

In order to obtain a result for the matrix element squared in the region S3 < A which is Lorentz- 
frame independent, the splitting function Pgg of the collinear cross section can be rewritten in terms 
of invariant variables according to 



1 



+ 



+ z{l-z) 



ko ■ h ko ■ kA 



+ 



+ 



k2-ki ' k2-k3 ' [p- (k-i + k4^)f 
u' te (54 + 55)2 



+ 0{9) 



(98) 



The expression ( p^ ) is Lorentz-invariant and can be used as an eikonal factor not only for the 
collinear configuration but also for the infrared region. Note that no additional divergences have 
been introduced in (|9^ ) since for /c2||A;4 and A;2||A;3 also the corresponding numerators vanish in the 
limit 9^0, k2\\k3 and A;2||A:4, respectively. Alternative ways of rewriting the splitting function in 
terms of invariants as e.g. 



1 



+ 0(9) 



(99) 



z p ■ ki 

only approximate the collinear matrix element squared but can not be used to describe the infrared 
region. 

An explicit distinction between collinear and infrared configurations as has been introduced in 
the centre of mass frame of the incoming partons according to (^TJ) is no longer necessary. Using 
the invariant formulation for the matrix element squared in the region S3 < A as given by the 
eikonal factor ( |98D , the phase space integrations necessary for the calculation of the corresponding 
cross section can conveniently be performed the centre of mass frame of the outgoing light partons: 
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,(100) 



By using ( pOD and (100) one finally obtains for the cross section of the process "yg ^ J/ip + gg in 
the region S3 < A: 
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(101) 



with Ce as defined in (| 

The cohinear matrix element squared for the process 'jg ^ J/i^ + QQ is found to be 
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with the gluon-(anti)quark sphtting function 
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By using the invariant expression 
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1 - e 



S5 



^1 +1- 



Si 



1 -e 



+ o{e) 



we find for the cross section of the process jg ^ J/'>p + QQ in the region S3 < A: 



(102) 
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(105) 



It has been checked explicitly that an expansion of the full matrix element squared in terms of 
small S3 leads to results identical to the equation above. 



C The kinematics of the photon-hadron cross section 

In this Appendix we discuss the kinematics of the hadronic differential and total cross section. The 
single-particle inclusive hadronic cross section for the process 

^{Ki) + P(Pi) ^ J/^(2p) + X (106) 

reads 

E /' . (107) 

dTidUl . ^ A^in X dtidui 

fg and /^--j denoting the parton distributions in the proton. The hadronic invariants are defined 
according to 

S = {Ki + Pif = s/x 
Ti = {Ki - 2pf - An? = h 

Ui = {Pi-2pf -Arri^ = ui/x . (108) 
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From (108) and the kinematical condition xS + Ti + xUi + 4m? > A one deduces 



A - Ti - 

It is convenient to rewrite the integration in terms of the J/tp energy variable z (|2|) and the 
J ftp transverse momentum p± . By using the relations 

s = xS ti = -- + Am^^ ui = -xzS (110) 



one finds 

dcT^^ 



E r dxfr{x,Q^)^-^ , (111) 



dp\dz 
where 

_ zA + pl + 4m^l-z) 

5^(1 -z) ■ ^ ' 

The total hadronic cross section is obtained by integration over pj_ and z: 

a-^{S,m^)= r^dz ff-dpl ^(S,pl,z,m^) , (113) 



mm 



with 

^max — 1 



Pimax = (l-^;)(5'z-4m2) -zA 

Pimin = . (114) 



Finally, the distribution is given by 



da^P{S,pl,m'^) _ /•^— dcjT-P 2 2 



dpi ^^^i„ dp'idz 

where 



d^ J 2 ^ {S,P±,z,m ) , (115) 



25 r V (5 + 4m2-A)2' 



mm 



The cut-off parameter A in the formulas above has to be included only for the calculation 
of the hard-gluon cross section. In our numerical program we have rewritten the In* A(i = 1,2) 
singularities of the virtual plus soft cross section as integrals over the momentum fractions of 
the partons x. The corresponding contributions have then been added to the hard cross section, 
cancelling the equivalent logarithms in this part so that the limit A — > could safely be carried 
out. 
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